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Summary 

Self- consistent solutions to the system of spinor and scalar field equations in General 
Relativity are studied for the case of Bianchi type-I space-time. It should be emphasized 
the absence of initial singularity for some types of solutions and also the isotropic mode 
of space-time expansion in some special cases. 



The aim of the paper is to find some exact self-consistent solutions to the equations, de- 
scribing spinor and scalar field system with the interaction Lagrangian Li nt = ip jCe ip ,a ^(S), 
<&{S) being arbitrary function of the invariant S = ipifj for Bianchi type-I space-time. 
Equations for partial choice of Q(S), while [$(S')]~ 1 = 1 + \S n , A being the coupling 
constant, n being some constant, have been thoroughly studied. It is shown that the 
equations, mentioned, can possess initially regular, as well as singular solutions, depend- 
ing on the sign of A, nevertheless singularity remains absent for solutions describing the 
field system with broken dominant energy condition. 

The Lagrangian for the interacting system of spinor, scalar and gravitation fields can 
be written as: 



R i 

L — h - 

2k 2 



- m ^ + 1^.^(5), (1) 



with R being the scalar curvature, k being the Einstein's gravitational constant. Function 
&(S) = 1 + \F(S), S = tfjipi describes the interaction between spinor and scalar fields, A 
being the interaction parameter. For A = the interaction vanishes and $(S) = 1. In 
this case we have the system of fields with minimal coupling. 
Bianchi type-I space-time metric can be chosen in the form [1] 

ds 2 = dt 2 - a 2 (t)dx 2 - b 2 (t)dy 2 - c 2 (t)dz 2 . (2) 

From Lagrangian (1) we will get Einstein equations, spinor and scalar field equations 
and components of their energy-momentum tensor. We will use Einstein equations for 
a(t),b(t) and c(t) in the form [1]: 



a a (b c 

- + - T + - 



a a\b c 



-«(?? - \t), (3) 

\t), (4) 



- + ? + - = -"fl?-H (6) 



a b c V u 2 

where point means differentiation with respect to t, and T = Tjf. 

Spinor and scalar field equations and components of its energy-momentum tensor can 
be written as follows: 

iy»V^ - + S)j, = 0, ®'(S) = ^, (7) 



(v^TyV^))=0, (8) 



1 d 
\f—g dx u 

T> = \g pu (fy M V„^ + $ lv V^ - V„fy^ - V„fy^) + W P *(S) - 6?L. (9) 
In (7) and (9) V M denotes covariant derivative of spinor, having the form [2]: 

= ^ - (10) 
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where T p (x) are spinor affine connection matrices. 7^(2;) matrices are defined for the 
metric (2) as follows. Using the equality 

9^(x) = e a ^{x)el{x)r] ah , ^(x) = e a ^{x)%, 

where n a b = diag(l, —1,-1,-1), % being flat space-time Dirac matrices, denoting a 
set of tetrad 4-vectors, we will get 

7° = 7°, 7 1 = 7V«W, 7 2 = 7 2 /^), 7 3 = 7 3 Mt). 
T^x) matrices are defined by the equality 

r,(x) = \g pa (x)(d,ey b -T^ s y i s , 

which gives 

r = o, i\ = \d(t)ff, r 2 = h(t)ff, r 3 = ^c(t)7 3 7°, (n) 

Flat space-time matrices we will choose in the form, given in [3]. 

We will study the space- independent solutions to spinor and scalar field equations (7), 
(8) so that 

V> = V(t), y? = <p(t). 
In this case solution to the equation (8) is: 

0W = ^y, C = const, r(t) = a(t)b(t)c(t). (12) 

In accordance with (12) the spinor field equation (7) can be written as: 

d f\ „ C 2 



if 



(- + i-)v-mV--P>(S)V = 0, (13) 



where P(S) = 1/$(S); P'(S) = % = -fj. For the components ip p = V p (t), p = 
1,2,3,4, from (13) one deduces the following system of equations: 

V r + ^-V r + i(m + ^-)v r = 0, r = l,2; (14) 



2r ' V 2r 2 
f „ ( C 2 P' 



Vi + 1 ^V l -i{m+ — 'jV l = 0, Z = 3,4. (15) 
From (14) and (15) we will find the equation for invariant function 

s = w = v*v 1 + v*v 2 - v*v 3 - v;v 4 ■. 



which leads to 



S+-S = 0, (16) 

T 

c 

S = — , Co = const. (17) 

T 
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As in the considered case P depends only on S, from (17) it follows that P(S) and 
P'(S) are functions of r = abc. Taking this fact into account, integration of the system 
of equations (14) and (15) leads to the expressions 



V r [t) = —=exp 



Vi(t) 



exp 



-i(mt + J Qdt) 
i(mt + J Qdt) 



r = l,2; 
I = 3,4. 



(18) 



where Q(t) 



C 2 P' 



~2~r~ i C r Ci - integration constants. 
Putting (18) into (9), we will get the following expressions for the components of the 



energy-momentum tensor for the interacting spinor and scalar fields 

-R, 



rpO 
J 



rp2 
1 2 



rp3 
J 3 



where 



N 



- 2 \cl + Cl-Cl-Cl){m + Q) 
r 



2iC n 



m + 



C 2 P 



2r 2 



R = 



C 2 



:l- 2 ( 



P + 



C P> 



i C 2 
T = T? = -N + — P - AR. 



Summation of Einstein equations (3), (4) and (5) leads to the equation 



T - = -k{tI + T 2 + T| - -r) = 3«( 



/mCo C C 2 P 



2t 



The first integral of the equation (22) takes the form: 

f 2 = 3K{mC T + -C 2 P + C 1 ), d 
2 

Final solution of the equation (22) reads 

dr 



At 3 



const. 



^/(mC r + \C 2 P + Gy 



±v3«(t + to); — const. 



(19) 

(20) 
(21) 

(22) 
(23) 
(24) 



Giving the explicit form of &(S), i.e. P — 1/$, from (24) one can find concrete function 
r(t) = a(t)b{t)c{t). Putting the obtained function in (18), one can get expressions for 
components of spinor function V p (t), where p = 1, 2, 3, 4. 

Let us express a, b, c through r. For this we notice that subtraction of Einstein equa- 
tions (3)-(4) leads to the equation 



a b ac be d fa b 
a b ac be dt \a b 



--T - + T + - =0. 



Equation (25) possesses the solution 



= Diexp 



(*/t)- 



D\ = const, X\ = const. 



(25) 



(26) 
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Subtracting equations (3)-(5) and (4)-(5) one finds the equations similar to (25), having 
solutions 

,14 , 1. . /■ ,14 . 

(27) 



a - = D 2 exp(x 2 h - = D 3 exp^X 3 J , 



where D 2 , D 3 , X 2 , X 3 are integration constants. There is a functional dependence between 
the constants D 1: D 2 , D 3 , X 1: X 2 , X 3 : 

D 2 = D ± D 3 , X 2 = X ± + X 3 . (28) 

Using the equations (26), (27) (28), we rewrite a(t),b(t),c(t) in the explicit form: 



a(t) = (DlD 3 )s T sexp 
b(t) = (D^D^rhxp 
c(t) = (D 1 D 2 



3 
X 1 



to 



X, 



to 



1 1 

3J 3 r 3 exp 



X l + 2X 3 



T~ l dt 



(29) 



to 



where to is the initial time. 

Thus the previous system of Einstein equations and interacting spinor and scalar 
field ones is completely integrated. In this process of integration only first three of the 
complete system of Einstein equations have been used. General solutions to these three 
second order equations have been obtained. The solutions contain six arbitrary constants: 
D\,D 3 , X\,X 3 and two others C\ and to, that were obtained while solving equation (22). 
Equation (6) is the consequence of first three of Einstein equations. To verify the cor- 
rectness of obtained solutions, it is necessary to put a(t),b(t) and c(t) in (6). It should 
lead either to identity or to some additional constraint between the constants. Putting 
a(t) , b(t) , c(t) from (29) in (6) one can get the following equality: 



1 

3r 



3f 



J 2 
2- + — 

T 6T 



X\ + X\X 3 + X 3 



-T 



(30) 



that guaranties the correctness of obtained solutions. 

To get the constant C 1 in (24) one can use the equation (30). Inserting f from (22), 
f 2 from (23) and 



J 



2 4 2 r 2 2r r 2 4r 3 



one deduces the identity, if 



c 1 = ±(x* + x 1 x 3 + xf), 



(31) 



which means that the constant C\ is a positive one. 

We will first study the solution to the system of field equations with minimal coupling 
when the direct interaction between the spinor and scalar fields remains absent, i.e. in 
the Lagrangian (1) ®(S) = 1. The reason to get the solution to the self-consistent 
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system of equations for the fields with minimal coupling is the necessity of comparing 
this solution with that for the system of equations for the interacting spinor, scalar and 
gravitational fields that permits to clarify the role of interaction terms in the evolution of 
the cosmological model in question. 

In this case the components of the energy-momentum tensor look: 

_ mCp C' , _ 2 _ 3 _ C' 
~ r + 2r 2 ' 1 2 3 ~ 2r 2 ' 



T = T: = ^-%, n + T] + Ti- 3 -T = - 3 -^. (32) 



Note that as the energy density Tq should be a quantity positively defined, the equation 
(32) leads to C > 0. The inequality C > will also be preserved for the system with 
direct interaction between the fields as in this case the correspondence principle should be 
fulfilled: for A = the field system with direct interaction turns into that with minimal 
coupling. 

Taking into account (32) equation (22) writes 

3 

f = -K,mC , (33) 

with the solution 

3 

r(t) = -KmC t 2 + Tit + t 2 , Ti,t 2 — const. (34) 

Putting r(t) from (34) into (18) and (29) one gets the explicit expressions for the 
components of spinor field functions V p (t) and metric functions a(t),b(t),c(t): 

V r (t) = ^e- imt , Vl (t) = ^e imt , (35) 



a(t) = {DlD^TzZ 1 
b{t) = (D^D^t^Z- 



c(t) = (D.Dly^Z-^, (36) 

where 

(t-hY 4 

\t-t 2 ) ' a 3«mC (ti - 1 2 ) ' 1 ' 



and ti j2 = — 3 ^ c ± 3K mCo V r i — 3/tmCor 2 are the roots of the quadratic polinomial in 
the right-hand side of (34). If the roots are real, i.e. if 

t\ - 3/tmC r 2 > 0, (38) 

the solution (34) is singular one, while in the opposite case it is not. Putting (34) into 
(30) one deduces the following relation between the constants: 

t\ - 3KmC T 2 = ^kC 2 + i [xl + XiX 3 + X^j . (39) 
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As the right-hand side of the equation (39) is positive, the quadratic trinomial in (34) 
possesses real roots and the solution obtained is singular one at initial time t — ti, whereas 
t\ > t 2 ,ti < t < oo. 

Let us study the solution (34)-(36) at t — > oo. Hence we have: r(t) ks jK,mC t 2 , 
and a(t) ~ b(t) ~ c(t) ~ t 2 ^ 3 , that leads to the conclusion about the asymptotical 
isotropization of the expansion process for the initially anisotropic Bianchi type-I space- 
time. 

Thus the solution to the self-consistent system of equations for the spinor, scalar and 
gravitational fields is the singular one at the initial time. In the initial state of evolution 
of the field system the expansion process of space-time is anisotropic, but at t — > oo there 
happens isotropization of the expansion process. 

To investigate the system of spinor and scalar field equations with direct interaction 
we will consider the partial case for choosing P(S): 

P(S) = 1 + XS n = 1 + A^|, (40) 

where A is the interaction parameter, n is some arbitrary constant. Inserting (40) into 
(24) one obtains 

dT = = y/zk, (41) 



mC T + \C 2 C%/2r n + C| 

where Cf = C 2 /2 + C\\ in (24) t has been taken zero, as it only gives the shift of the 
initial time. 

Let us study different cases of choosing A and n. I. A > 0, n > 0. In this case (41) 
leads to the following behavior of r(t): 

3 

at t — > oo r(t) ks -KinCot 2 — > oo, (42) 



at t^O r(t) 




2 



x/2+l 

-> 0. (43) 



Note that (42) coincides with (34) at t — > oo. It leads to the fact that in the case 
considered, the asymptotical isotropization of the expansion process of initially anisotropic 
Bianchi type-I space-time takes place without the influence of scalar field. In this case 
the initial state is singular: r(0) = 0. 

Thus, the evolution of the interacting fields system at A > and n > is qualitatively 
the same as that of the system with minimal coupling. 

II. A = —a 2 < 0, n > 0. The equation (41) takes the form: 



/ 



dr 



^mC T - a 2 C 2 C%/2T n + C\ 



VSKt. (44) 



From (44) follows: 



3 

at t — > oo r(t) ~ -nmCot 2 — > oo, 



i.e. as well as in the previous case the asymptotical isotropization of the expansion process 
of initially anisotropic Bianchi type-I space-time takes place. But r = cannot be reached 
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as in this case the denominator of the integrand in (44) becomes imaginary at r 
There exists the minimum value r m i n = t > 0, which is defined from the equation 



mCor^ 1 + C 2 2 r« 



0. 



It means that for A < and n > there exist regular solutions to the previous system of 
equations. The absence of the initial singularity in the considered cosmological solution 
appears to be consistent with the violation for A < 0, of the dominant energy condition 
in the Hawking- Penrose theorem [1]. 

III. A > 0, n — —k 2 < 0. In this case the equation (41) takes the form: 



dr 



mC T + \C 2 r k2 /2Cf + Cl 



3nt. 



(45) 



Let us study concrete solutions for some values of k 2 . 
a) k 2 = 1. Then from (45) one gets: 



M = m + 

2C 2 



AC 2 
2C 2 



(46) 



The solution (46) is singular one at initial time t = and asymptotically isotropic, 

b) k 2 = 2. The equation (45) writes 



dr 



sjmC Q r + \C 2 T 2 j2Cl + a 



3nt. 



(47) 



Integration of (47) leads to 

r(t) 



r 2 
AC 2 



'V3k\C 
, V2C 



t — mCn 



(48) 



A = ^2\C 2 C 2 /Cl - m 2 Cl 
From (48) one gets: r(t ) = 0, where to is defined from the equation: 



Ash\ 



V3k\C 
. V2C 



t - mC = 0, 



(49) 



i.e. the solution (48) is singular at initial time t = t . 
At t -> oo 



(50) 



The solution (48) describes initially (i.e. at to) singular and asymptotically (i.e. at t — > oo) 
isotropic Bianchi Type-I cosmological model. Note that in this case the transition to the 
isotropic regime happens exponentially. 



IV. A = —a < 0, 



n 



—k < 0. In this case the equation (41) takes the form: 
dr 



^mCoT - a 2 C 2 r k2 /2Cf + C. 



3nt. 



(51) 



Let us consider concrete solutions for some values of k 2 as in III. 
a) k 2 = \. In this case one gets from (51): 



/ \ a 2 C 2 
Vr + y/nln | yfr - yjr{ | ) = VSrtt, y/n = th, (52) 



which leads to 

at t -> oo r(f) « — -> oo, (53) 



( ° 2C ' 2 \ 

at t -> -oo r(f) -> n = ^^3/2 ) • ' : ' ! 1 ; 



From (53) and (54) one comes to the conclusion that the solution (52) is initially (i.e. at 
to = — oo) regular one and at t — > oo asymptotically isotropic, 
b) k 2 = 1. In this case (51) writes: 

di~ i — . s 

VSnt. (55) 



(m-a 2 C 2 /2C^jC r + C. 



If in (55) m — a 2 C 2 /2Cq > 0, then the solution coincides with (46), where M = m — 
a 2 C 2 /2C 2 . At m - a 2 C 2 /2C 2 = —T 2 < from (55) one gets 



In this case r(t) possesses 



r(t) = ^ - 3 -^t 2 . (56) 



maximum at t = 0, i.e. r(0) — 



T 2 C n ' 



2C 2 

and minimum at t\ 2 = T — — z.e. r (^i 2) = 0. (57) 

v3/tT 2 C 

The solution obtained describes the cosmological model, which begins to expand at t±, 
acquires its maximum at t — and then collapses into a point at t 2 . 
c) k 2 = 2. From (45) one gets: 

/ dT = y/zk. (58) 

J JmC T - a 2 C 2 T 2 /2C$ + Cf 



(59) 



Integrating (58), for r(t) one gets the following expression: 

aCy3nt^ 



1/2 



mC + Asin 



V2C , 

where A = ( m 2 Cl + 2a ^ C2 ] . From (59) follows that r(£ ) = 0, where 



r = ^= — arcsm[ — — I, (60) 



then acquires maximum 

r{t m ax) = ^[mC + AJ, 

where 



tmax 2V^aC' 
and further at t — ti again turns to zero: r(ti) = 0, where 

V2C . (mCo 

t\ = 7r H == arcstnl — - — 

V^aC V A 

Thus the solution (59) describes the cosmological model, which begins to expand at 
to, acquires its maximum at t max and then collapses into a point at t±. 
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